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Abstract. 

We study the properties of short gamma-ray bursts (GRBs), based on the 
assumption that they are all connected to the binary neutron star mergers, whose 
formation mechanism contains a large amount of uncertainty. In particular, the merger 
time distribution is focused on, which is essential for short GRB rate estimate; we 
simply parameterize it as Pm{t) oc with lower cutoff time-scale r. The short GRB 
rate, calculated using such a parameterization with a specific model of cosmic star 
formation rate, is found to be quite insensitive to the value of r, but considerably 
sensitive to a. With such formation rate models, the peak flux distribution observed 
by BATSE is fitted. In the fitting process, the local short GRB rate po is treated as 
one of free parameters, giving the best fit when po = 0.1-2 yr~^ Gpc~'^, depending on 
the other free parameters as well as the formation rate model represented by (t, a). 
We also discuss several implications for the intrinsic properties of possible sources of 
short GRBs, such as jet structure, intrinsic luminosity function, and intrinsic duration 
distribution. 
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1. Introduction 

In recent years, our knowledge of gamma-ray bursts (GRBs) enormously improved 
through many observations of GRBs themselves (using gamma-rays) as well as their 
transient component, i.e., afterglow, which can be detected in various wavebands from 
radio to X-rays. Especially from the afterglow observations, their cosmological origin 
has been estabhshed, which revealed that GRBs are very luminous astrophysical objects 
whose isotropic energy emitted by gamma-rays ranges 10^^-10^^ erg. Furthermore, 
evidence of jetted GRBs arises from long term radio observations and from 
observations of achromatic breaks in the afterglow light curves |2]. Using the observed 
break time, the jet opening angle and intrinsic gamma-ray energy were inferred 0, and 
the authors have concluded that the gamma-ray energies, after correcting the beaming 
effect, are clustered around ~ 10^^ erg. However, the remarkable progress is restricted 
only to the long duration class of GRBs, whose time duration Tgo exceeds 2 s, where Tgo 
is the interval over which 5-95% of the burst counts accumulate. On the other hand, no 
afterglows have been observed thus far for short duration GRBs (Tgo < 2 s), and as its 
consequence, there is no direct knowledge of their redshifts and intrinsic luminosities. 

Because of its energy and time scales, GRB has been considered to be connected to 
massive compact stars, such as supernova explosions and binary neutron star mergers. In 
fact, several observations indicate that long GRBs are likely associated with supernova 
explosions the best evidence for the connection is the spectroscopic discovery of 
SN 2003dh in the afterglow spectrum of GRB 030329 jSj- The absence of afterglow 
observation for the short GRBs, on the other hand, restricts our discussions concerning 
their origin. However, it is natural to interpret that the phenomenon that triggers short 
GRBs is different from that triggering long GRBs; this is due to clear bimodality of the 
distribution of Tgo accumulated in the BATSE data 0. In the present paper, we assume 
that the short GRBs are associated with the binary neutron star mergers, i.e., neutron 
star-neutron star (NS-NS) or neutron star-black hole (NS-BH) coalescence. Numerical 
studies actually suggest that short duration burst can be explained by the NS-NS or 
NS-BH merger models \!\. 

Coalescence rates of binary compact object system have been estimated in two ways: 
(1) theoretically, based on the predictions of binary population synthesis calculations 
(e.g., [Hlinilin]) or (2) empirically, based on the observed sample of galactic binary pulsars 
(e.g., [H]). At present, both methods appear burdened with significant uncertainties: 
the theoretical approach due to the many poorly constrained evolutionary model 
parameters and the empirical estimates due to small number sample of observed NS- 
NS systems. Thus, both theoretically and observationally estimated coalescence rates 
considerably change as a result of discovery of (theoretical) new merger formation 
channels or (observational) tight binary systems; in fact, there were such discoveries 
in both approaches. In the theoretical population synthetic approach, a new formation 
channel was found to be relevant in reference p2j, resulting in the increase of the number 
of tight binary compact object systems that possibly merge in less than 1 Myr from 
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their formation. From the observational aspect, a highly relativistic binary neutron star 
system has been recently discovered |TH], which will merge in about 85 Myr as a result 
of gravitational wave emission, and the observational rate estimate increased by a factor 
of 6-7 compared with the previous one |14j . 

As we assume in the present study, if the short GRBs are truly associated with 
binary neutron star mergers, it is expected that the observations of short GRBs and 
their afterglows would provide rich information on the picture of binary systems, such 
as formation rate history and surrounding environment. In this paper, we investigate 
the formation rate history of short GRBs using the peak flux distribution obtained by 
the BATSE detector, particularly focusing on the merging time distribution (for an 
earlier approach, see reference It depends on the physical processes of the binary 

formation that are still quite uncertain as already discussed in the previous paragraph. 
Several past studies [ini HZl EE] have also used the peak flux distribution, but in order 
to obtain implications for cosmic star formation history or luminosity functions (LFs) of 
long/short GRBs. In the present paper, we rather aim at deriving physical properties of 
NS-NS or NS-BH mergers from the short GRB observation, assuming the connections 
between them and that the observationally inferred star formation rate (SFR) correctly 
represents nature. Although there remain huge amount of uncertainties concerning the 
cosmic SFR history, many observations with various wavebands seem to be roughly 
consistent with each other at low-redshift region after some correction law of dust 
extinction is applied. 

Because the available data are quite restricted, owing to absence of any afterglow 
signals, we believe that it would be better to use simple parameterization for merger 
time distribution as described in section |2l It is found that the GRB rate history is 
quite sensitive to the merger time distribution. In section El we briefly model necessary 
properties of short GRBs such as the LF and photon spectrum. The available BATSE 
data as well as its efficiency are also summarized in the same section. Section|3]is devoted 
to showing the results of parameter fitting. The best fit values for each parameter as 
well as its allowed region at la level are presented. Finally in section IHl we discuss the 
obtained properties concerning neutron star mergers as a potential source of short GRBs, 
by comparing our result with those obtained by other approaches such as population 
synthesis; we discuss the intrinsic GRB properties such as lumnosity and jet structure. 
In addition, the redshift distribution of BATSE short bursts is inferred using the best 
fit parameters, and then we derive the intrinsic duration distribution by correcting time 
dilation effect. 

2. Formation history of short gamma-ray bursts 

By recent progressive observations at various wavebands, the cosmic SFR is fairly well 
known, although there remains debate on especially high-redshift universe. As our 
reference model, we use a simple functional form for the SFR per comoving volume, 
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which is given in reference |18] (labeled as 'SF2'), as 

exp(3.42) + zf + VLkY'^ ^ ^ 1.. s 

and we adopt the standard ACDM cosmology {Q^ = 0.3, = 0.7, Hq = 70 h^Q km s~^ 
Mpc~^); figure n shows SFR-z relation with this model (solid curve). The SFR sharply 
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Figure 1. Formation rate liistory of siiort GRBs. Eacli rate iiistory is obtained with 
the assumption that the short GRBs are associated with binary neutron star mergers. 
Dot-dashed, dashed and dotted curves correspond to different lower cutoff time-scale 
of coalescence distribution Pm{t) with its slope fixed to be a = —1.0. The absolute 
values are normalized to best fit the observed peak flux distribution as discussed in 
section^] Cosmic SFR density, from which GRB rates are generated, is also shown as 
solid curve. 



increases from the current epoch z = to 2; ~ 1.5, and then it keeps the same level 
to farther past universe, z = 5. Behaviours at high-redshift region are quite uncertain 
because the galaxy LF has not been established yet there as well as the amount of dust 
extinction is highly ambiguous. Therefore, it may be useful to compare several SFR 
models which differ at high redshift from each other, as actually done in reference JH] 
(for long GRBs). In the present study, we further adopt two other SFR models also 
given in reference ^H] as 'SFl' and 'SF3', and compare the results obtained by these 
models to that by the standard model SF2. The models SFl and SF3 are given as a 
functional form by 

exp(3.4z) \n^(l + z)^ + QaY^^ ^ ^ 1 . . 
exp(3.82;) -|- 45 (1 -|- zy^'^ 
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RsYz{z) = 0.22/170 



exp(3.05z - 0.4) + zf + Vt^f'^ 

exp(2.93^) + 15 (1 + 2)3/2 



Mq yr"^ Mpc 



-3 



(3) 



Behaviours of these two functions are roughly the same as that of SF2 at low-redshift 
region z < 1.5, but on the other hand, they significantly differ at 2; > 1.5; at that 
region, the SFl and SF3 models give exponentially decreasing and increasing formation 
rate, respectively, compared with the constant behaviour of SF2. In fact, as shown in 
section|l[ it is found that the results obtained by the SFl and SF3 models are essentially 
unchanged compared with that by SF2, since the bulk data should be contributed rather 
low redshift region z < 1.5, at which three SFR models are only slightly different. 

The short GRB rate, RcRBit), can be calculated assuming its association with 
binary neutron star mergers, by using the formula 



where tp represents the formation epoch of galaxies, about which we assume z{tp) = 5. 
Pm(t) is the probability distribution function of merging time of the binary system 
from its formation. This distribution can be calculated theoretically using population 
synthesis method [El El CHI, t>ut resulted in various different distributions depending on 
studies, which means that quite large ambiguities should be accompanied. In the present 
paper, therefore, we take more simplified approach; i.e., a simple parameterization such 
as Pmit) oc with lower cutoff time-scale r is adopted. This is because we beheve that 
this approach would be more useful to obtain some intuitive implications. We adopt 
as the power law index a values (—0.5, —0.7, —1.0, —1.5), and (0.1, 1, 20) as TMyr, where 
TMyr represents the lower cutoff time-scale in unit of Myr. One of the parameter sets 
we adopt, i.e., (a,rMyr) = (—1.0,20), was also investigated in reference [12]. As already 
mentioned in section Q a new formation channel was pointed out in reference ^21, with 
which a large amount of tight binaries, which merge within 1 Myr from their formation, 
possibly form. It resulted in decrease of lower cutoff time-scale r; for instance, merging 
time distribution calculated in reference pQ] is found to peak at ~ 0.3 Myr, which is 
significantly shorter than the previously considered value, i.e., > 10-100 Myr. 

Figures ^ and |21 shows the short GRB rate density as a function of redshift due 
to binary neutron star mergers that is calculated with equation (j3]). Normalization of 
the absolute value is obtained by the parameter fitting that is discussed in section HJ 
Dependence on the parameter r is shown in figure^with a = —1.0. As the value of r is 
increased, the local rate evolution to z ~ 1 becomes shallower than that of original SFR 
model, from which the short GRB rate was generated via equation This is because 
the time delay effect from the formation to coalescence is more prominent in the case 
of larger r. On the other hand, in figure |21 we compare models with different values for 
power law index a assuming TMyr = 1; the GRB rate history is found to change more 
dramatically depending on the value of a. Using smaller a, which increases the relative 
contribution of tight binaries to the total rate, the GRB rate becomes steeper at lower 
redshift; especially for a < —1.5, the GRB rate history becomes quite good tracer of 
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Figure 2. The same as figure ^ but evaluated for various values of a with fixed r to 
be 1 Myr. 

cosmic SFR. In consequence, the short GRB rate history is quite sensitively dependent 
on the unknown parameters of binary neutron star systems. 

3. Peak flux distribution and theoretical models 

We define short GRBs as those with Tgo < 2 s, and take from the BATSE current burst 
catalogue4 The BATSE peak fiux is given in three different time-scales of 64, 256 and 
1024 ms. We use the 64 ms triggering timescale because the average duration of short 
GRBs is ~ 0.3 s, and we have to use a shorter timescale than this for the real peak fiux. 

The peak fiux distribution can be theoretically calculated using the following 
formulation: 

which is in unit of sr^^ yr~^ cm^ s. Here, dV/dz is the comoving volume element per 
unit solid angle, ip{L) is the peak LP of short GRBs, e(P) is the detector efficiency as a 
function of photon fiux, and the factor (1 + accounts for cosmological time dilation. 
In our notation, the burst peak luminosity L is given as a function of redshift z (that 

X http://cossc.gsfc.nasa.gov/batse/ 



Short GRBs and neutron star mergers 



7 



corresponds to distance) as well as the observed peak flux P. The comoving volume 
element is simply represented by the following formula with cosmological parameters: 



dV^ 
dl 



where di^{z) is the luminosity distance to an object at redshift z. The efficiency e(P) 
is for 64 ms trigger timescale taken from the BATSE current burst catalogue. In the 
catalogue the trigger efficiency is given only above P ~ 0.8 cm~^ s~^, though the data 
exist at lower flux range. Hence, we restrict data fltting only to this high flux region, 
where the trigger efficiency is given in the BATSE catalogue. 

The relation between peak luminosity L and observed peak flux P (and also with 
redshift z) is obtained if we assume the differential rest-frame photon luminosity of the 
source S{E) (in units of s~^ keV~^). As S{E), we use single cutoff power-law model 

i/3 + 2)E 



S{E) oc E^ exp 



peak 



(7) 



with observationally inferred value for short GRBs |19j, i.e., (3 = — 0.58,£'pcak = 355 
keV. The peak photon flux P observed at the Earth in the energy band -Emin < E < -Emax 
and emitted by an isotropically radiating source at redshift z can be written by 

„ (i+^)ji£^g::g(g)dg 
^ Sdjw • 

where we adopt -Emin = 50 keV and -Emax = 300 keV, which are adjusted to the energy 
band of the BATSE detector. On the other hand, an 'isotropic equivalent' burst peak 
luminosity is deflned in the photon energy 30-2000 keV as 



2000 keV 



ES{E)dE, 



(9) 



30kcV 



from which the relation connecting luminosity with peak flux and redshift, L{P, z), can 
be derived, combined with equations ((Tj) and (jHl). 

We simply assume that the burst luminosity distribution does not evolve with 
redshift and adopt simple functional form for ip{L), 



L 



exp 



L 



(10) 



Lor(-7 - 1) 

where L denotes the isotropic peak luminosity in the 30-2000 keV energy range at burst 
rest frame, as deflned by equation ©, 7 is the asymptotic slope at the bright end, Lq 
marks a characteristic cutoff scale. We label this peak LF as single power-law (SPL) 
model. For comparison, we also adopt another LF whose functional form is represented 
by 
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Table 1. Models used for parameter fittings. 
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where we assume Ai = 30 and A2 = 10. In contrast with the previous SPL-LF 
[equation (fTUj) ]. this LF is labeled as double power-law (DPL) model. Further as 
two power-law indexes 61 and ^2, we choose two sets of their values, i.e., (^1,^2) = 
(—0.1, —2.0), (—0.6, —3.0), and name each model as DPLl and DPL2, respectively; this 
choice of the power-law index was actually adopted in the recent publication [201, which 
focused on the long GRBs, on the contrary to the present study. For both LFs, a proper 
normalization 'ip{L)dL = 1 is ensured. All the information given above allows us to 
calculate theoretically the peak flux distribution using equation 0, and to compare it 
with the observational data points. 

In consequence, we prepare 28 models, for which the parameter fitting is performed; 
each model is characterized by (1) the adopted values for merger time distribution a and 
r, (2) the adopted LFs and (3) cosmic SFRs. We summarize all the models considered 
from this point on in table ^ and we label each model by the model number 1-28 as 
shown in the table. 
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4. Results of parameter fitting 

For each model given in tabled we set the best fit value of several free parameters, i.e., 
Pq,'J,Lq in the case of SPL-LF and po,L^ in the case of DPL-LF, where po represents 
the local {z = 0) short GRB rate. Observational data selected as short GRBs using 
the criterion Tgo < 2 s are converted into rates per unit time per unit solid angle by 
considering the effective coverage of ~ 0.48. We show in figure El the peak flux rate 
distribution with the best fit theoretical models calculated by equation For figures 
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Figure 3. Differential short GRB number counts as a function of peak photon flux. 
Curves represent the best fit models obtained from the parameter fitting to the data 
points. In each panel, several models (definitions are given in table ^ are shown: (a) 
models 3, 7 and 11, (b) models 5-8, (c) models 7, 15 and 19. Panels (d)-(f) are the 
same as (a)-(c) but the fit is restricted to the 12 data points from the highest flux 
value [compared with 15 points in the panels (a)-(c); see text for more details]. 
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Table 2. Best fit parameters and one sigma errors for models 1-12. They are 
characterized by the different values of (r, a). As for other properties, SF2 and SPL-LF 
are adopted. The values of the minimized as well as goodness of fit are shown in 
the last two columns; the degree of freedom is 12. 
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ini^a)-(c), data points were used for fitting when their flux is larger than ~ 0.8 cm^^ s^^, 
above which the trigger efficiency is given in the BATSE current catalogue; there, 15 
data points exist and thus it results in degree of freedom (z/) of 12 and 13, in the case 
of SPL-LF and DPL-LF, respectively. 

FigureEfa) shows the result of parameter fitting for the models 3, 7 and 11; i.e., we 
investigated the dependence on the lower cutoff time-scale by fixing a to be —1.0. The 
dependence on the power-law index of merger time distribution a, on the other hand, 
is shown in figure Efb) with fixed value of TMyr = 1, i-e., comparison among models 
5-8. The best-fit values to three free parameters, pcT and Lq, as well as their allowed 
ranges at 68% confidence level are summarized in table |21 for models 1-12. In the same 
table, we also show the values of as well as a goodness of fit (g.o.f.) in the last two 
columns. The value of is not sufficiently small compared to the degree of freedom, 
but the fit is still acceptable. As clearly seen from table the dependence on the lower 
cutoff time r is very weak, while the values of the best-fit parameters are rather sensitve 
to the power-law index a. These characteristics obviously reflects the strong (weak) 
dependence of the short GRB rate on a (r), which is already discussed in section |2l and 
shown in figures [T] and |21 Table El suggests that at present, with any models we can 
obtain reasonable fit to the data; no models are highly favoured or disfavoured compared 
to the others. 

Several implications for the burst parameters, i.e., po)7 and Lq, are also obtained 
by the fit. The local short GRB rate is found to be 0.13-1.3 yr~^ Gpc~^, which 
is comparable to or in most cases larger than the local long GRB rate inferred by 
the analysis in reference ^Hl; i.e., 0.12 yr^^ Gpc^^, by some factor to one order of 
magnitude. Although only 1/3 fraction of total number of detected GRBs is short type. 
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Table 3. Best fit parameters and one sigma errors for models f3-20. They are 
cliaracterized by the adopted LFs and different values of a. As for other properties, 
SF2 and TMyr = 1 are adopted. The values of the minimized as well as goodness of 
fit are shown in the last two columns; the degree of freedom is 13. 
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Table 4. Best fit parameters and one sigma errors for models 21-28. They are 
characterized by the adopted SFRs and different values of a. As for other properties, 
SPL-LF and TMyr = 1 are adopted. The values of the minimized as well as goodness 
of fit are shown in the last two columns; the degree of freedom is 12. 



Po Lq 
Model SFR a (yr^^ Gpc^^) ^ i^iqSI gj.g g-i) ^2/^ g ^ £ 



21 


SFl 


-0.5 


1 9 +0.3 
^■^ -0.2 


-6 2+^-2 


4 ci+3-O 
^•^-2.4 


17.0/12 


0.16 


22 


SFl 


-0.7 


0.72±°:- 


-3 

"3-^-1.2 


9 0+3.6 
'^•°-2.3 


17.3/12 


0.14 


23 


SFl 


-1.0 




_2 q+O-8 


1 q+30 


17.4/12 


0.13 


24 


SFl 


-1.5 


^■-•^^-o.oe 


-2.6;gj 




17.4/12 


0.13 


25 


SF3 


-0.5 


1.3 ±0.3 


-6.4t- 


4.6t|i 


16.7/12 


0.16 


26 


SF3 


-0.7 


n 77+0.26 
' ' -0.16 


-4 2+^-2 


q q+2.9 
0.0_2 1 


17.0/12 


0.15 


27 


SF3 


-1.0 


^■■J^-o.ie 


q 1+0.6 


2.5ti8g8 


17.2/12 


0.14 


28 


SF3 


-1.5 


^■^"-0.17 


-2.8t;:^ 


9 7+10.5 
^- ' -2.3 


17.3/12 


0.14 



the current estimated formation rate of short GRBs dominates over that of long GRBs. 
The remaining parameters related to the LF can compensate this apparent conflict. The 
best fit values for such parameters are —6.4 < 7 < —2.7, and (2.1-4.6) x 10^^ erg s~^ for 
Lq. According to reference ^Hj, on the other hand for long GRBs, the best fit values 
for these parameters are Lq = 6 x 10^^ erg s~^ and 7 = —2.9, with which the luminosity 
distribution of long GRBs is expected to be brighter and flatter (owing to larger Lq 
and 7) than that of short GRBs. Thus, smaller values for Lq and 7, obtained for short 
GRBs, are expected to push more GRBs below the detection threshold, reducing the 
total number of detected GRB rate. Therefore, the resulting number of short GRBs 
is smaller than that of long ones, although the local rate itself po is dominated by the 
short population. 

We further show the same results but for models 13-20 and 21-28 in tables El and 
EJ respectively. The best-fit flux distribution for the DPL-LF models (model 15 and 
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Table 5. The same as table El but the fitting is restricted to the rather high-flux 
region. The degree of freedom is reduced to 9. 

Po Lo 



Model 




Ol 




/ 


flO^i erg s~i) 


1 V 


f 


1 


20 


-0.5 


1.5 ±0.4 


7 9+2.7 
' -^-3.2 




8.9/9 


0.44 


2 


20 


-0.7 


1.2 ±0.4 


-6. 8+33 


4.8;3i 


8.5/9 


0.49 


3 


20 


-1.0 


^•"3J^-0.19 


-3.8t°:^ 


9 q+3.1 


8.9/9 


0.45 


4 


20 


-1.5 




_o Q+0.6 
'-'•"-40.2 


9 q+173.2 
^■"-•-2.0 


9.0/9 


0.44 


5 


1 


-0.5 


1.4 ±0.4 


7 0+2.7 
' -^-3.2 




8.9/9 


0.45 


6 


1 


-0.7 


1.1 ±0.4 




4.813-1 


8.4/9 


0.49 


7 


1 


-1.0 


0.38t°i? 


-3.4±^8^4 


9 R+2i0.4 
^■"-2.6 


9.0/9 


0.44 


8 


1 


-1.5 






9 q+41.8 

^■'-'-2.0 


9.0/9 


0.43 


9 


0.1 


-0.5 


1.4 ±0.4 


7 9+2.7 
' -^-3.2 




8.9/9 


0.45 


10 


0.1 


-0.7 


1.1 ±0.4 






8.5/9 


0.49 


11 


0.1 


-1.0 


n QQ+0.32 
U.OO_o.24 


q 9+0.7 
•-'•^-47.9 


9 r+207.5 


9.0/9 


0.43 


12 


0.1 


-1.5 


16+019 

'^■^"-0.12 


_9 q+0.5 


9 q+40.8 
^■"-•-2.0 


9.0/9 


0.43 



19) are compared with that for the SPL-LF model 7 in figure Efc). In the case of the 
DPL-LF models, the quality of the fit becomes slightly better, but not to the extent that 
we can conclude that the short GRB data strongly favours DPL luminosity distribution. 
The adopted SFR does not give any significant difference from one another; this clearly 
indicates that at present the local contribution z <l.b from SFRs dominates the bulk of 
data. Because many observations with various wavebands give roughly consistent results 
with the SFl, SF2 and SF3 models at low-redshift region, i.e., the SFR increases from 
the current epoch 2; = at least to z = 1, we believe that the uncertainties concerning 
the adopted SFR are small at low redshift and therefore giving no significant ambiguity 
to the parameter values as found from the comparison between tables El and |31 

Until this point, we used all the available data given in the BATSE current burst 
catalogue. However, it should be noted that because the trigger efficiency used in the 
present paper does not include the correction concerning the effects of atmospheric 
scattering. As the result, the efficiencies at low fiuxes are underestimated as also 
mentioned in the BATSE catalogue itself, although we do not know whether our 
calculations given above are affected by such an effect. In order to investigate that 
point, we repeated the same procedure given above, but with restricting the fit to the 
12 data points from the high values of peak fiux (compared with the 15 data points 
previously adopted), and its result is summarized in figures intd)-(f), tables and IHl 
This time, the degree of freedom is reduced to 9 for SPL-LF and 10 for DPL-LFs. It 
is noticeable that the quality of the fit becomes considerably better for both SPL and 
DPL models. We believe that this fact suggests that the fitting improvement is not due 
to the wrong modeling (at least of LFs) but rather due to exclusion of lower-fiux data 
points, which may be affected by the atmospheric scatterings. Although the values of 
the best-fit parameters are quite consistent with the previously obtained values in tables 
El and El the accompanying errors become significantly large. Further, every model gives 
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Table 6. The same as table 13 but the fitting is restricted to the rather high- flux 
region. The degree of freedom is reduced to 10. 



Model LF a (yr"i Gpc^^) (lo^i erg s'^) xV^ g-O-f- 



13 


DPLl 


-0.5 


1 S +0-3 
^■^ -0.2 


2 0+°-3 


10.9/10 


0.37 


14 


DPLl 


-0.7 


1.3 ±0.2 


2.1 ±0.3 


8.9/10 


0.54 


15 


DPLl 


-1.0 


A2+°-^° 

^•^^-0.10 


3.91^-? 


8.7/10 


0.56 


16 


DPLl 


-1.5 




6.2t^-? 


8.9/10 


0.54 


17 


DPL2 


-0.5 


^•L* -0.4 


q 7+0.8 

-J-'-i.g 


9.8/10 


0.46 


18 


DPL2 


-0.7 


-■^•D -0.4 


q 7+1.3 
-J-'-l.? 


8.3/10 


0.60 


19 


DPL2 


-1.0 


n 07+0.09 
^■•J' -0.12 


9.3±?° 


9.6/10 


0.48 


20 


DPL2 


-1.5 






8.6/10 


0.57 



reasonable fit at almost the same level; no models can be excluded or restricted by the 
present analysis. 

5. Discussion 

5.1. Implications for jet structure and intrinsic luminosity 

It is of interest to compare these properties with those inferred from other 
theoretical/observational approaches. The local merger rate of binary NS-NS or NS-BH 
sysyem calculated by population synthesis code ranges by many orders of magnitude. 
For instance according to reference P], it ranges 0.01-80 Myr~^ galaxy"^ depending 
on the various parameters, with 1 Myr~^ galaxy"^ in the case of standard choice of 
parameter set. On the other hand, more recent study ^0] claims that larger number of 
NS-NS merger can be created; the coalescence rate ranges 1-300 Myr^^ galaxy^^ with 
standard value of 50 Myr~^ galaxy^^ 

In this section, we compare these merger rates obtained from their standard models, 
i.e., 1 Myr~^ galaxy"^ [9 or 50 Myr~^ galaxy"^ ^U] with our inference of po using short 
GRB observations. We simply convert the rate per galaxy into that per volume assuming 
that the number density of galaxies is ~ 10-2 Mpc-^ Using this conversion factor, the 
rate 1 (50) Myr"^ galaxy"^ corresponds to 10 (500) yr~^ Gpc~^, with which the values 
Po we have obtained in the previous section should be compared. From tables EHEl we 
take Po = 0.5 yr^^ Gpc~'^ as our reference value, which is significantly smaller than 
those derived from the population synthesis approach; our evaluated value is 0.05 and 
10~^ of the local merger rate from references jH] and PU]; respectively. However, this 
disagreement may be a natural consequence of the jet structure of GRBs or of the fact 
that not every merger produces GRBs. Assuming that all the binary mergers containing 
at least one neutron star certainly produce short GRBs, we can suggest the beaming 
angle of the jet using the rate difference obtained above. Considering that the fraction 
of apparent GRB rate is only 6'^/2 of the intrinsic GRB rate (coalescence rate), the 
derived values for the jet opening angle 6 are 18° and 2.6°, which are derived from the 
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Table 7. Comparison with other theoretical/observational rate estimates. The 
galactic merger rate is converted into global rate with an assumption that the local 
galaxy density is 10"^ Mpc"'^. As the local short GRB rate pq and cutoff luminosity 
Lq, the values given in table El are adopted with the specific model 7. 



Method 


Reference 


Merger rate 
(Myr-i Gal-i) 


Jet angle 

en 


r int 
^0 

(lO'*^ erg s-i) 


Population synthesis 


Fryer et al 9 


1.2 


13 


5.7 




Belczynski et al JOI 


52.7 


2.0 


0.13 


Observation 


Kalogera et al Jl] 


180 


1.1 


0.038 



ratio of po obtained from our approach to that from the local coalescence rate taken 
from references P| and [lOj, respectively. These inferences of the jet opening angles are 
quite well consistent with those for long GRBs determined by the achromatic break of 
the afterglow light curves 3j. The intrinsic values corresponding to the luminosities, 
Lq and L^, are suggested to be also smaller than the values summarized in tables |2Hni 
by a factor of 0.05 or 10~^, depending on whether the opening angle is 18° or 2.6°. 
We summarize these particulars in table Q but in which the specific model 7, i.e., 
(riviyr, a) = (1, —1-0), SF2 and SPL-LF, is used for the evaluation with the values given 
in table |21 In the same table, we also show the values obtained by comparing our po 
with the observationally inferred local coalescence rate, which was recently found to be 
rather large as 180]'^^44 Myr^^ galaxy"^, although the estimate contains a fair amount 
of uncertainty jjlj. In this case, the obtained jet angle is a little bit small but is still 
marginally consistent with that of long GRBs observed by the achromatic break of 
afterglow light curves. In consequence, although there is no direct evidence supporting 
the connection between short GRBs and binary neutron star mergers, we can explain 
observed data quite naturally with such a connection as well as jet-like nature of short 
GRBs. 

In the above discussion, we have assumed that the jet structure is uniform in the 
opening angle, i.e., the kinetic energy does not depend on the solid angle. On the other 
hand, several models that assume structured iet, i.e., the kinetic energy depends on the 
solid angle, are proposed (e.g., see references [211 122]), and whether the jet is uniform or 
structured is matter of controversy [2II1 ES] • At the end of this subsection, therefore, we 
consider the possible case of structured jet model; we adopt the power-law distribution 
of kinetic energy per unit solid angle, i.e., e(^) oc and k = 2 as the canonical 
value [21]. A remarkable property of the structured jet model is that the isotropic LP 
is determined to be ip{L) oc L^^^'^^'^, because the apparent luminosity depends on the 
direction alone j2Ill221- In our calculation, this corresponds to the case of SPL-LP model 
summarized in table |21 Prom the table, the canonical value of 7 = —2, expected from 
the structured jet model with = 2, is not favoured, while the statistical significance 
is not large. Steeper jet structure (due to larger k) gives flatter LP (larger 7) and it is 
further more disfavoured. In addition, our results are also against the quasi- universal 
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Gaussian jet structured model, which predicts ip{L) oc L ^ even if some dispersion of 
parameter values is allowed [221 121] ■ 

5.2. Inferred redshift distribution and intrinsic burst duration 

Although we do not have any redshift information for short GRBs, we can infer their 
redshift distribution detected by BATSE, by using the formulation. 

Figures m and El show the redshift distribution dN/dz calculated using equation (|T^ 
with the GRB formation rate Rgkb{z) shown in figures ^ and |21 In figure |3J we also 




Figure 4. Redshift distribution of short GRBs detected by BATSE. Three models 3, 
7 and 11, which are different by the value of lower cutoff time-scale r with a fixed to 
be —1.0, are compared with each other. The distribution of long GRBs obtained by 
parameters given in reference |18 | is also shown for comparison. 



show the expected redshift distribution of long GRBs detected by BATSE using the 
best fit parameters taken from reference ^H], for comparison. From the figure, the short 
GRBs are considered to be more localized relative to the long ones, but its degree is 
relaxed when the slope of the merger time distribution becomes steeper (for smaller a). 
These redshift distributions would be another probe of its origin if the accompanied 
afterglow data were accumulated. 
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Figure 5. The same as figure01 but evaluated for models 5-8, i.e., with various values 
of a with fixed r to be 1 Myr. 



Although we classified the bursts with Tgo < 2 s as short GRBs, the observed 
duration is affected by cosmological time dilation. In order to obtain intrinsic duration 
distribution for short GRBs, we use an assumption that it obeys log-normal distribution 
as well as the redshift distribution of the short GRBs already obtained via equation (fT^ . 
The rate of detected bursts per unit redshift range dz per unit intrinsic duration dAtjnt 
is obtained by 

d^N = P(Atint)dAtint^dz , (13) 

where P(Atint) represents the probability distribution function of intrinsic duration, 
which we assume log-normal type; dN/dz has already been obtained in the previous 
paragraph. This relation can be rewritten as a function of observed duration At = 
(1 + z)Atint and redshift z, and then be integrated over redshift as: 



dAt .In dz dAt 

' P{{1 + z)-^At) dN 



dz 

. dz, (14) 
1 + z dz 

with which we can obtain the values for central duration as well as the standard deviation 
of the log-normal distribution P(Atint), by fitting the observed duration distribution by 
BATSE. As the result we found that (Atint) = 0.24-0.31 s and aiogAw = 0.50-0.55, 
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depending on the adopted parameters of (r, a). 
6. Conclusions 

Recent progresses concerning GRBs such as evidence for colhmated jet structure or 
association with supernova explosions are all restricted to long duration GRBs with 
Tgo > 2 s. For the short GRBs, on the other hand, because of the lack of any afterglow 
information, we only poorly understand them. However, a clear bimodal duration 
distribution and several numerical simulations suggest that the origin of short GRBs is 
different from that of long ones, and the binary neutron star merger event remains as a 
viable candidate. 

Based on the assumption that all the short GRBs with Tgo < 2 s are connected 
to the binary neutron star mergers, we analyzed the current available data taken from 
the BATSE catalogue. In the analyses, we parameterized the merger time distribution, 
which is essential for short GRB rate estimate, as Pm(^) with lower cutoff time- 
scale r. The short GRB rate, calculated using such a parameterization with equation 
(0)) and generated from the model for cosmic star formation rate SF2 in reference |18j . 
was found to be quite insensitive to the value of r, but significantly sensitive to a, as 
shown in figures ^ and El With such formation rate models, the peak flux distribution 
observed by BATSE was fitted. We further adopted two different functional form for 
isotropic luminosity distribution, i.e., SPL and DPL. In the fitting process, we used 
several free parameters, po? 7? -^o (Po; -^*) in the case of SPL-LF (DPL-LF) model, where 
Po represents the local short GRB rate density and the others are the parameters related 
to the LFs. The result of parameter fitting is summarized in figures 0] as well as tables |21- 
m We found that every model gives acceptable fit with reasonable parameter values, and 
that they are essentially insensitive to the adopted SFR models (SF1-SF3) because the 
bulk of the data points are expected to be attributed to local redshift region z < 1.5. The 
resulting values for the local GRB rate are found to be 0.1-2 yr~^ Gpc~^. Furthermore, 
the quality of the fit considerably improves when we omit several data points at low peak 
fiux as shown in tables El and (HI This may be attributed to the fact that the BATSE 
triggering efficiency is underestimated at such a region, because it does not contain the 
correction to the atmospheric scatterings. 

We also discussed several implications for the intrinsic properties of possible sources 
of short GRBs. (1) By comparing with the galactic coalescence rate of binary NS-NS 
systems obtained by theoretical population synthesis or observations of binary neutron 
star systems, and also by assuming that all such mergers certainly produce short GRBs, 
we constrained the jet structure as well as the intrinsic luminosity function. The 
obtained values for a specific model is summarized in table (2) The structured jet 
model, where the kinetic energy depends on the direction, e.g., e(^) oc 9~'', was also 
considered; a remarkable characteristic of this model is that the isotropic LF becomes 
SPL function 'ipiL) oc L'^'"^^^. According to our analyses, the model with > 2 is 
disfavoured although the statistical significance is still small. (3) Although the redshift 
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distribution of short GRBs is not known yet since there is no discovery of afterglow 
signals, we can infer it using equation (fT^ and best fit parameters, as shown in figures 
m and El With this distribution and the observed duration distribution of short GRBs, 
we derived the parameters for intrinsic duration distribution which we assumed to have 
a log-normal property, by correcting the effect of cosmological time dilation. 
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